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Abstract 

The magnetic properties of the ferromagnetic Kondo-lattice model 
(FKLM) are investigated. Starting from an analysis of the magnon spec- 
trum in the spin-wave regime, we examine the ferromagnetic stability as a 
function of the occupation of the conduction band n and the strength J of 
the coupling between the localised moments and the conduction electrons. 
From the properties of the spin-wave stiffness D the ferromagnetic phase 
at zero temperature is derived. Using an approximate formula the critical 
temperature Tc is calculated as a function of J and n. 



1 Introduction 

The intensively investigated Kondo-lattice model (KLM) describes the exchange 
coupling of the itinerant electrons to a system of permanent, localised magnetic 
moments. These moments have their origin in the partially filled 5/-, 4/- and 
3rf-shells of materials such as actinides, rare earth elements, transition metals 
and their compounds. 

The electrons of the conduction band mediate an indirect exchange interaction 
between the local moments, which order collectively below a critical temperature 
Tc- The Hamiltonian of the Kondo-lattice model is defined as follows llll: 



ija i 

The first term describes the hopping of the conduction electrons between dif- 
ferent lattice sites at and Rj. Ty are the hopping matrix elements, which 
represent the amplitude of the respective hopping process. The hopping integral 
is 



^^^■-^E<k)e^>^(^-^^), (2) 

k 
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where e(k) is the Bloch dispersion. The operators af^ and aja denote the cre- 
ation and annihilation operators for conduction electrons with spin a. 
The second term describes the exchange coupling between the system of local 
moments and the conduction electrons. This term is similar to the Heisenberg 
Hamiltonian |^ and describes the local character of the interaction, ai is the 
spin operator of the conduction electrons, and Si the operator of the localised 
spins. 

A very important fact is the sign of the coupling constant J. If J < 0, the itin- 
erant electrons and the local spins couple antiferromagnetically. One result is 
the existence of a nonmagnetic phase caused by screening of the local moments 
[|[ The term "Kondo-lattice model" implies a negative J. If J > the 
two subsystems couple ferromagnetically. Because the model with positive J 
differs from the KLM only in the sign of J it is refered to as the "ferromagnetic 
Kondo-lattice model (FKLM)" or ''sf-(sd-) model". 

Interest in the FKLM has been on the increase following the rediscovery of 
'colossal magnetoresistance' in various manganites, where a large change in re- 
sistivity associated with the FM transition occurs. The first observation of 
this effect has been made in 1950 Manganites have a perovskite structure 
Ri-xXa;Mn03, where R=La, Pr, Nd and X=Sr, Ca, Br, Pb. The manganite 
Lai_j;Caa;Mn03 is a prototype for the double exchange model 0|, which can 
qualitatively describe ferromagnetism in manganites. However, it cannot repro- 
duce the complex phase diagram of Lai-^CaajMnOs, where antiferromagnetic, 
ferromagnetic and paramagnetic phases as well as phase separation has been 
observed 

A procedure to reproduce this phase diagram using the FKLM has been pre- 
sented in a paper by Dagotto et a/.jsj These authors calculate phase diagrams 
for one- and two-dimensional systems with different total spins S. Both the 
existence of different magnetic phases and phase separation were demonstrated. 
A similar phase diagram has been calculated for infinite dimensions ||^. It is 
generally believed, therefore, that the FKLM is a good starting point for the 
determination of a phase diagram of manganites even in three dimensions. This 
is plausible if the t2g electrons of the Mn^+ ions in Lai_a;Ca2,Mn03 are assumed 
to be quasilocalised S = |-spins and the eg electrons are assumed to be itin- 
erant electrons moving through the S = ^ lattice and being ferromagnetically 
coupled to the localised moments. 

Another class of substances which are described by the FKLM are the ferro- 
magnetic metals Gd, Tb, Dy and doped EuX ||l^. Gd is an excellent prototype 
since its total spin of | is the largest possible for the rare earth elements. The 
main focus of this paper are the ferromagnetic metals. In the framework of the 
FKLM the main difference between these materials and the manganites is the 
value of of the couplings strength J, which is one order of magnitude larger in 
manganites. The coupling constant in ferromagnetic metals is in the order of 
0.1 eV |l|. 

The main aim of this work is to determine the ferromagnetic phase in three 
dimensions at T = K by analysing the stability of the spin-wave spectra. 
The magnon energies will be calculated by determining the poles of a magnon 
Green's function. The same approach has been used by Wang Q, whose nu- 
merical results showed two significant properties of spin-wave spectra in the 
FKLM, namely the vanishing of the acoustic mode in the continuum and the 
effect of an anomalous softening of the spin-wave dispersion at the edge of the 
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Brillouin zone. The second effect tias been experimentaffy demonstrated in tlie 
manganite Pro.esSro.syMnOs |Q. The author of did not, however, investi- 
gate systematically the spin- wave dependence on the coupling strength and the 
band occupation n. The dependence of the spin wave energy on the coupling 
strength in the FKLM has been calculated by Furukawa for a fixed band oc- 
cupation. An increase of the magnon energy with increasing coupling strength 
has been shown ||l3|] . 

In this work a systematic investigation of the spin wave energy at different band 
occupations and coupling strengths is presented. From the spin wave stiffness D 
we can determine the ferromagnetic phase in a J — n diagram in 3D at T = K. 
It will be of particular interest to see whether and at what parameter configura- 
tion anomalous softening occurs. Finally we calculate the critical temperature 
dependence on J and n using an approximate formula taken from the Heisenberg 
model. 



2 Theory 

In this section we derive the magnon Green's function for low temperatures in 
the random phase approximation (RPA). We then discuss briefly the resulting 
energetic structure of the magnon system. 
We use the identities for the electron spin operators (h—l) 



and for the local moment spin operators 

sf = list + s-y. sf = ^ist S-) (4) 

to find a more useful form for the s/-term in Eq. (j^) 

^'^Z = -^^i^'^S■n^a + S^a+_^a^a) (5) 

i.{7 

with a =]= +, 0" =[= — , Z| = 1 and Z| = —1. The S'f are the spin flip operators 
which cause spin deviation in the local moment system. 

The first term on the right-hand side of Eq.(H) is an " Ising-term" , since it 
describes the interaction of the z-components of localised spins Q . The second 
term is non-diagonal in the spin indices and is called the "spin-flip term" . 
We want to investigate the properties of our system at very low temperatures. 
In this limit we can apply the spin-wave approximation using the Holstein- 
Primakoff transformation ]l6| . This expression can then be simplified for low 
temperatures. By applying this approximation to one gets: 
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(6) 

bf ,bi are Bose operators. The elementary excitation of the spin system can 
be derived from the retarded magnon Green's function {{S^ ; S^))e, which is 
defined as foUows: 



{{S+; Sr))E = -i / dt expi^iEt) {[S+{t), ^"(O)]-). (7) 







(...) means thermodynamic averaging, while [...]- denotes the commutator. 
In the spin wave approximation, the magnon Green's function is given by 

{{S+;S-))E = 2S{{bf,b+))E. (8) 

The equation of motion of ((6,^; b'^))E is 



Em-M))E = + 6+)) 



E 



On the right-hand side of the equation of motion we have two higher Green's 
functions. The spin-flip Green's function represents a characteristic part of 
the interaction between electrons and magnons since a magnon is created or 
annihilated by an electron spin-flip process. These processes are of primary 
importance in determining the magnetic properties of the FKLM. We therefore 
decouple the Ising-function only: 



{{b^lM,■b-^^))E = (10) 

The spin-flip function is treated more carefully. We deflne a general spin-flip 
Green's function Tik,j{E): 

T,kAE)^{{at^aki;b+))E. (11) 

If we now apply Eq.(|lO|) to the equation of motion (^ and use the Fourier 
transform Gq{E) — ((5q;6+))£; of Gij{E) ~ {{bi\b1j))E we get the following 
form of the equation of motion in q-space: 



EG^{E) = l + J{a^)G^{E)^^V2SY,e-^<^^-^^^T,,AE) (12) 

= ^E^-("k.). (13) 
k.cr 
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We now derive the equation of motion for Tikj{E): 



ET.kAE) = ([a,;tafex,fe+]_) + (([a,tafci,iJ]_;6+))B 

= '^^{Tk-nX in,j{E) — TniTnkjiE)) + JSTikj{E) 
n 

- ^ (^{{72ka+.^aki;b+))E + {{nia+aki;b+))E^ (14) 

- ^V2S [{{bka+akr,b+))E - {{b.a+aki;b+))E) . 

On the right-hand side of Eq.(p^ we have a set of higher Green's function 
which describe more complex interaction processes between itinerant electrons 
and magnons. These higher Green's functions are RPA decoupled, preserving 
spin and particle conservation, as foUows: 



((nfea+aA;|;6+))£; w {nk)Tik,j{E) 

{{bka+akr,bp)E ~ {a+ak^)GkjiE). (15) 

Applying the Fourier transform of the spin-flip Green's function to ^ we and 
using Eq.(|^) we obtain the following approximate expression for the equation of 
motion: 



r ((nkT)-(nk+qu))Gq'(i?) . . 

''^'''^ ' 2x/iV i?-e(k + q')+6(k)-J(5-^Ep(V)' 

At this stage we can obtain an analytical expression of the niagnon Green's 
function by solving Eq. (p^ : 



- E-J{.^)-J^Sxici,E) ^''^ 

where 



^i^' ) 2Ar^i?-e(k + q) + e(k)-J(^-iEp(-P» ^ ^ 

and (rip) is the occupation number of the magnons. At very low temperatures 
the number of magnons in the system is very small since the system is almost 
saturated ferromagnetically. We therefore set (rip) = for numerical simplicity. 
In ||l^ the magnon Green's function ((S'+(— q); 5^ (q))) is derived by solving 
the equation of motion for finite temperatures in RPA. If we set (5^) = S in 
this more general expression, our result with (rip) — and q = — q is obtained. 
In the RPA theory, the magnon energy is renormalised by the term xi^.iE), 
which is due to electron-magnon interaction. This term is given explicitly by 
the band structure, the occupation of the electron system and the coupling 
strength. 
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In this work we wish to investigate the energetic structure of the magnons. The 
magnitude of the dispersion is influenced by the real part of the magnon self 
energy, while the magnon lifetime is described by its imaginary part. The self 
energy Mq{E) can be defined by the formal solution of the retarded magnon 
Green's function 



The magnon self energy in RPA can now be derived from Eg. (p!?)) : 

Mq{E) - J(a^) + J^S{Rexiq, E) + z/mx(q, E)). (20) 
By applying the Dirac identity we can calculate the real and imaginary part of 

x(q,^): 

/mx(q, E) = ^ E(("kT) - ("k+qi))^(£; - e(k + q) + e(k) - JS) (22) 

k 

where V denotes the principal value. This gives rise to two collective excitation 
modes and a one particle continuum as discussed in pTf . 

Outside this continuum the principal value at the right-hand side of eq. (^ij) is 
equal to the function itself, the imaginary part vanishes and the magnon lifetime 
becomes infinite. The poles of the Green's function, i.e. the collective excitation 
energies, are then defined by the implicit equation 



J(ct^) - j25i?ex(q,S) = 0. (23) 

The stability of the local moment fcrromagnet at low temperature is deter- 
mined by the collective excitations from the ferromagnetic ground state. The 
excitations into the acoustic mode are certainly dominant. We therefore limit 
ourselves to an analysis of the dependence of the acoustic mode on J and n. 
In order to solve Eq.(p3|) an expression for the (riko-) is required. It is very com- 
mon to treat the electronic system in a mean-field approximation (see e. g. p^). 
This leads to a situation in which the electron-magnon interaction is included 
in the investigation of the magnon system, but not of the electronic system. In 
this work we wish to investigate how the ferromagnetic properties of our model 
system are determined by the RPA theory. We therefore limit ourselves to a 
mean- field treatment of the electronic part. In this approximation the expecta- 



tion values of the k-dependent occupation numbers in Eq. (21) is 



^"''"^ " e/3(.(k)-M-i J.„s) + 1 ■ (2^) 
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In this work we assume the system to be a simple cubic lattice with nearest- 
neighbor interaction only. The electronic dispersion of such a system is given 

by 

d 

e(k) = -2t^cos(fc,) + io, h e[-7r;7r], t > 0. (25) 

i=l 

The midpoint of the band to is chosen to be zero. The relation between the 
band width W and the hopping term t is Vl^ = 12t. 

3 Results and Discussion 

In this section we discuss the properties of the spin- wave dispersion Eq.(p3|) 
with increasing band occupation n at different values of the coupling J. The 
coupling of the conduction electrons to the local moments leads to an indirect 
interaction between the spins. 

In Fig. |l| spin wave spectra in the [100] direction are sketched for J = 0.2eV 
and J = 0.25eV at different band occupations n. By fiUing the conduction 
band with electrons the magnon energy increases up to a maximum at a given 
band occupation n. This maximum value increases with increasing J. After 
reaching the maximum value close to the quarter-filled band, the magnon en- 
ergy decreases monotonically with further increasing band occupation. 
The excitation energy of a spin-wave mode is the energy needed to cause a de- 
viation of the spin system of unity from the ferromagnetically saturated state. 
Hence the larger the spin-wave energy the larger the ferromagnetic stability. 
This is of course a qualitative statement which must be substantiated by quan- 
titative calculation of, for example, the critical temperature Te- 
la Fig. ^ we have sketched the behaviour of the magnon energies by further 
increasing the band occupation. The spin-wave spectra become negative. 
The existence of negative excitation energies is of course unphysical since the 
magnons are excitations from the ground state. There is thus no stable ferro- 
magnetic ground state beyond a critical band occupation n^. The determination 
of Uc requires the definition of a stability condition. The system becomes ferro- 
magnetically unstable if the magnon energy becomes negative. 
A widely used parameter for the ferromagnetic stability is the spin-wave stiff- 
ness D. In it has been shown that the expression for the magnon energy 
can be rewritten as -^(q) = Dq^ for small q. When D becomes zero the sys- 
tem becomes ferromagnetically unstable. In the intermediate coupling regime 
0.2eV < J < 0.7eV, this is indeed a unique criterion since the spin-wave energy 
becomes zero at small q. 

However, an interesting effect can be observed at very small J. In Fig. || the 
magnon spectra for a system with J — 0.1 eV are sketched. At n = 0.24 the 
spin-wave dispersion starts to decrease significantly for large q while the spin- 
wave stiffness remains almost the same. This effect is known as 'anomalous 
softening'. A second minimum appears in the dispersion curve at the righthand 
edge of the Brillouin zone. At n = 0.255 the spin-wave energy at the X-point 
becomes zero. Hence at small J a ferromagnetic instability exists at large q. 
This effect has been observed experimentally in the manganite Pro.eaSro.syMnOa 
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Q. But in manganite systems a strong coupling between the conduction elec- 
trons and the localised moments exists. It is thus not very clear, whether or not 
the anomalous softening found in our calculation is comparable to the experi- 
mental results. A more complex theory which assumes some rf-character of the 
conduction electrons has already been presented in order to explain the mea- 
sured anomalous softening ||l8|. A reasonably good agreement between theory 
and experimental data has been found. Despite this unsolved problem anoma- 
lous softening is certainly a model property of the FKLM as well. A very 
significant consequence of anomalous softening is the fact that the spin-wave 
stiffness D is almost constant while the magnon energy at large q is lowered. 
The decrease in magnon energy is equivalent to a decrease in ferromagnetic sta- 
bility and hence of the critical temperature Tc- In our theory these facts would 
lead to an increase of the relation S- . This effect has also been observed experi- 
mentally in manganite systems [19|. In this paper a ferromagnetic ground state 
is assumed. It is thus only possible to report a unique ferromagnetic instability 
at J = 0.1 eV. The cause of the anomalous softening at low J and hence the 
nature of the new phase is currently the subject of further research. 
We have defined above the critical band occupation ric where the spin-wave 
stiffness becomes zero and the system ferromagnetically unstable. By calculat- 
ing Tic at different coupling J we can sketch the ferromagnetic phase in a J — n 
diagram, shown in Fig. |4[ 

The ferromagnetic phase occupies a large part of the phase diagram. The criti- 
cal band occupation increases by increasing the coupling strength J. There is no 
ferromagnctism around half filling (n — 1). One gets a qualitative explanation 
of this result by looking at the energy of different spin states of the conduction 
electrons at half filling. All lattice points are occupied with one electron. If 
all electrons are in the spin-up state, no virtual hopping between the lattice 
points is possible. Virtual hopping reduces the energy of the system. Thus the 
stable spin configuration corresponds to the state with the maximal hopping 
amplitude, i.e. where all nearest neighbour electrons have opposite spin. Since 
we assume a ferromagnetic coupling between the two electron subsystems, such 
a spin configuration corresponds to an antiferromagnetic order of the local mo- 
ments. 

At J = 0.1 eV we have sketched two points. The black point represents the 
band occupation where the magnon energy vanishes at the right-hand edge of 
the Brillouin zone. The second point is sketched at ric- Hence the procedure 
described to determine the phase border is no longer unique. 
The structure of the ferromagnetic phase of our system is very similar to those 
calculated in |^ and These calculations were performed using a variety of 
different methods, such as quantum Monte Carlo calculations, for systems with 
different dimensions d and total spins S. It is therefore quite problematic to 
compare these phase diagrams directly. However the similarity between these 
results is remarkable. The qualitative behaviour of the ferromagnetic phase 
seems to be a general property of the FKLM. 

Finally we have calculated the critical temperature using an approximate for- 
mula which can be derived from the spin- wave approximation of the Heisenberg 
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model (see e. g. 0): 

where E{q) is the spin- wave dispersion of a local moment system with direct 
exchange interaction. 

This formula extrapolates the critical temperature from the magnon energy of 
the system in the spin-wave regime. A quantitatively better calculation of the 
critical temperature would require a self-consistent calculation. 
At this stage we can only give a qualitative description of the dependence of 
the critical temperature on the band occupation n and the coupling strength 
J. In our system the exchange interaction is mediated by the conduction elec- 
trons. The FKLM-Hamiltonian can be mapped onto an effective Heisenberg- 
Hamiltonian pT[ | . Information about the kind of interaction is then buried in the 
effective exchange integrals J(q), which are functionals of the electronic self en- 
ergy. This description leads to the modified RKKY-interaction Following 
the same idea we can describe the excitation of the local spins in the spin- wave 
approximation as a system of noninteracting magnons: 



i/ = i?o + Ef^eT/(q)^qV (27) 

q 

E(q)=ri^^"^(q) is the spin-wave dispersion (p3|). The renormalised spin-wave 
energies are equivalent to the effective exchange integrals in the modified RKKY- 
interaction . In our theory we can therefore rewrite the approximate formula 
(Uw for the critical temperature. 




This equation is used to calculate the dependence of the critical temperature by 
evaluating the magnon energies in the whole Brillouin zone for given n and J. 
Plots of Tc versus n for different J are sketched in Fig. ||. From Eq. ( ^ ) it is 
evident that the functional dependence of Tc is equivalent to the n-dependence 
of the magnon energy. The critical temperature increases by increasing the 
band occupation up to a maximum. After passing this maximal value the criti- 
cal temperature decreases monotonically before it vanishes at the critical band 
occupation ric (black points). Since the error of the numerical evaluation of Tc 
becomes very large near the critical band occupation, the value of nc was taken 
from the phase diagram in Fig. ^ The good agreement between the Tc versus 
n graphs and the respective Uc values supports the stability condition using the 
spin-wave stiffness D defined above. Good qualitatitive agreement with the Tc 
dependence on n reported in Jill has been found. This agreement is perhaps 
unexpected, since a more realistic description of the electronic system is used 
in This could be due to the fact that the exchange integrals used in Q 
are to first order very similar to the expression of the real part of x(Q5 ^) i^i 
Eq.(pl|) which renormalises the expression of the magnon energy Eq. (p3|) with 
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respect to electron-magnon interaction. The qualitative behaviour of the criti- 
cal temperature is mainly determined by these interaction processes. It is thus 
reasonable that the inclusion of those processes in the perturbative analysis of 
the magnon system leads to a good description of the qualitative behaviour of 

An exception is the Tc{n) dependence at J = O.leV. The calculation stopped 
at n = 0.21. Above this band occupation we have to sum over two singularites 
in Eq.(|2^) due to the effect of anomalous softening. A numerical calculation is 
now rather complicated and was not done at this stage. However we sketched 
the critical band occupation Uc = 0.38, where the spin-wave dispersion at small 
values of q vanishes. If we assume the same qualitative dependence of Tc as 
found for the larger J this ric is too large. Hence the stability condition using 
D does not work at this J. 

A similar qualitative dependence of the critical temperature, on the band oc- 
cupation, including a maximum at quarter filling, has also been found for the 
strong coupling limit (J — > oo) of the FKLM assuming classical spin (5 = 1) 
and using Monte-Carlo methods pO|] . 

Another interesting result concerns the critical temperature as a function of the 
coupling strength J in Fig. ^. It is sketched for three different band occupations 
n. The Tc versus J dependences for n = 0.1 and n = 0.06 are quite similar. The 
critical temperature increases with increasing coupling strength J. This result is 
intuitively correct. For n = 0.4 we found a critical coupling strength Jc- Below 
Jc, Tc is zero and the system is not ferromagnetically ordered. The existence 
of a Jc has also been shown in ||ll| for higher band occupations. A different 
qualitative behaviour compared to ||ll| has been found at large J. While in 
that work J remains constant above a certain value, we did not find saturation. 
This difference could be caused by the different electronic self energies used. We 
have not calculated the critical temperature for J < O.leV. For this region we 
expect the normal RKKY-behaviour Tc ^ J^, as shown in pT[ |. 



4 Summary and Conclusions 

We have investigated the ferromagnetic properties of the FKLM through an 
analysis of the spin-wave spectrum at very low temperature. The spin-wave 
approximation was applied to the model Hamiltonian and the magnon Green's 
function was derived within the RPA approximation. The numerical analysis of 
the magnon energy spectrum was limited to the acoustic branch of the collective 
excitations. 

The results of n- and J-dependent calculations of the magnon spectra have been 
presented. The magnon energy for fixed coupling strength J increases with in- 
creasing occupation of the conduction band up to a maximum value close to 
quarter filling. After reaching the maximum value the energy decreases as n 
increases and eventually vanishes. A stability condition for the ferromagnetic 
state was defined by using the spin- wave stiffness D. If the magnon energy at 
small q vanishes at the critical band occupation Uc, i.e. D becomes zero, the 
system becomes ferromagnetically unstable. The ferromagnetic phase in three 
dimensions has been determined by calculating Uc for different J. There is 
a good agreement to the ferromagnetic phase obtained for FKLM-model sys- 
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terns with different total spin S and dimension d in and g]. The critical 
band occupation Uc increases by increasing the coupling strength. There is no 
ferromagnetism around n = 1. The good agreement between the different calcu- 
lations may lead to the conclusion that the structure of the ferromagnetic phase 
is a general property of the FKLM. 

In addition anomalous softening is found at J = O.leV. The method to deter- 
mine ric is thus no longer unique. It is not yet possible to link this result to 
the anomalous softening observed experimentally. This problem as well as the 
physical reason of this phenomenon is an interesting subject of future research. 
Finally the critical temperature dependence on J and n was presented. A qual- 
itative agreement to the results shown in ||l^ was found. A more careful treat- 
ment of the electronic system and a selfconsistent calculation of the critical 
temperature is currently the subject of further research. 
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Figure 1: Spin- wave dispersion E{q) in [100] direction for a system with S = 3.5, 
band width W = 1.0 eV at T = K. The dispersions for two different coupling 
constants J = 0.2 eV and J = 0.25 eV at different band occupations n are 
shown. On the left-hand side of the figure the magnon energy increases with 
increasing n, whilst at the right-hand side it decreases. 
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Figure 2: Spin- wave dispersion E{q) in [100] direction for a system with S = 3.5, 
band width H^ = 1.0eVatT = 0K. The dispersions for two different coupHng 
constants J = 0.2 eV and J = 0.25 eV at different band occupations n are 

shown. 
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Figure 3: Spin-wave dispersion E{q) in [100] direction for a system with S = 3.5, 
band width W = 1.0eVatT = 0K. The dispersion for different n at J = 0.1 
eV are shown. 



0.7 
0.6 
0.5 
^0.4 
0.3 
0.2 
0.1 





0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

n 



Figure 4: Ferromagnetic phase at T = K, 5 = 3.5, and W = 1.0 eV 
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Figure 5: Critical Temperature Tc as a function of n at different coupling con- 
stants J. 
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Figure 6: Critical Temperature Tc as a function of J at different band occupa- 
tions n = 0.06, 0.1, 0.4. 
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